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O ■ ABSTRACT 
O 
(N 

Qh' We report the detection of redshifted iron Ka absorption lines in the Chandra 

ly^ • LETG spectrum of the narrow-hne quasar, PG 1211+143. The absorption hues 

are observed at 4.22 keV and 4.93 keV in the quasar spectrum, corresponding to 
4.56 keV and 5.33 keV in the rest frame of PG 1211+143. From Monte Carlo 
simulations, the chance probability of both lines being false detections is low 

■ at 1.36 X 10~^. Highly redshifted ionized iron Ka (Fe xxv or Fe xxvi) is the 

■ most plausible identification for the lines at their observed energies. If identified 
Q"^ ■ with H-like iron Ka at 6.97 keV, then the relativistic velocity shifts required are 
^ ■ 0.40c and 0.26c. The extreme velocities can be explained by pure gravitational 
O ■ redshift if the matter exists in a stable orbit within 6 gravitational radii of the 
^ ■ black hole. This would require a Kerr metric for the black hole. Alternatively 

the absorption may be the result of matter infalling directly onto the black hole, 

O 

eg ■ This matter may originate in a failed outflow or jet, which does not escape the 

gravitational potential of the black hole. 
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1. Introduction 

One prime motivation of high energy research into Active Galactic Nuclei is to find 
evidence for the putative massive black hole, through the gravitational redshift imparted 
on photons emerging within a few gravitational radii of the event horizon. To date, most 
of the effort has focused on the iron Ka emission fine, predicted to show a relativistically 
broadened profile (Fabian et al. 1989; Laor 1991), which appeared to be common among 
Seyfert 1 galaxies, as observed by ASCA (Tanaka et al. 1995; Nandra et al. 1997). However 
the situation emerging from the newer XMM-Newton and Chandra data is more ambigu- 
ous, the spectra showing that narrower (non-relativistic) lines are commonplace (Yaqoob & 
Padmanabhan 2004). Indeed the red-wing of the broad iron line is often hard to distinguish 
from the effects of the warm absorber (Reeves et al. 2004; Turner et al. 2005a). 

Given the ambiguity in modeling the broad iron line, additional diagnostics of the inner 
disk material are needed. Narrow, transient and redshifted iron emission lines of modest 
statistical significance have recently been claimed in several AGN (Turner et al. 2002; Dovciak 
et al. 2004; Porquet et al. 2004) possibly originating from local hot-spots on the disk surface. 
Potentially these lines can be tracked to reveal orbital motions in the inner disk (Dovciak 
et al. 2004; Turner et al. 2005b). However httle evidence has been found to date for infall 
of matter onto the black hole. In this paper, we present evidence of highly redshifted iron 
Ka absorption lines in the narrow-line Seyfert 1, PG 1211-1-143 with Chandra LETG. This 
may either be interpreted as infall onto the massive black hole, or gravitational redshift 
occurring within < 6Rg of a Kerr black hole. PG 1211+143 (at z=0.0809) is a prototypical 
'Big Blue Bump' quasar, and one of the brightest AGN in the soft X-ray band (Walter 
& Fink 1993). A previous XMM-Newton observation of PG 1211+143 showed blue-shifted 
absorption from a high velocity {v ~ 0.1c) outfiow (Pounds et al. 2003). The soft X-ray 
properties of PG 1211+143, which confirm the blue-shifted absorber, will be presented in a 
subsequent paper. Values of Ho=70kms~^ Mpc~^, qo—0-0 and Aq = 0.73 are assumed and 
errors are quoted at 90% confidence (Ax^ = 2.7), for 1 parameter of interest. 



2. Chandra and XMM-Newton Observations of PG 1211+143 

PG 1211+143 was observed by Chandra with the Low Energy Transmission Grating 
(LETG) in 3 consecutive orbits between 21 June 2004 and 26 June 2004 of approximately 
45 ks exposure per orbit (see Table 1). An XMM-Newton observation was conducted si- 
multaneously with the first Chandra orbit, but due to high and variable background, this 
was reduced to only ~20ks useful exposure. The ACIS-S detector was used on the focal 
plane, with the S3 chip being at the telescope aimpoint. Chandra LETG source spectra. 
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background spectra and response matrices were generated using CIAO 3.2.1. Streaks were 
removed from the ACIS-S CCD array, while events with grades 0, 2, 3, 4, and 6 were se- 
lected. After screening an exposure of 133.6 ks was obtained. Spectra and responses from 
the ±1 orders were extracted and combined to maximize signal to noise. The time-averaged 
combined first order LETG count rate was 0.31 counts s~^. Data from the zeroth order 
image were not used, due to photon pile-up. The 1st order spectra were binned so that one 
bin represents the FWHM resolution of the LETG (AA = 0.05), ensuring sufficient counts 
per bin to use the minimization spectral fitting technique. 



3. LETG Spectral Analysis 

To parameterize the X-ray continuum, the LETG spectrum was fitted with a broken 
power-law model, modified by neutral Galactic absorption of iVn = 2.8 x lO^*' cm~^ (Dickey 
& Lockman 1990), over the 0.4-8.0 keV bandpass. The soft and hard X-ray photon indices, 
break energy and model normalization were allowed to vary; the best fit values are Fsoft = 
2.77 ±0.03, Lhard = 1.96 ±0.05, break energy E = 1.7 ±0.1 keV and normalized flux of 1.6 ± 
0.1 X 10~^ photons keV~^ cm~^ s~^ (at 1 keV). The integrated flux is 1.42 x 10~^^ erg cm~^ s~^ 
(0.4-8 keV), corresponding to an unabsorbed X-ray luminosity of 2x10^^ ergs~^. The spectral 
fit obtained is formally unacceptable, with x^/dof = 686/544 (dof=degrees of freedom) for 
a rejection probability of 2.9 x 10^^. This is due to several ionized emission and absorption 
lines (e.g. from O, Ne) below 1 keV, although the continuum is fitted well above 2 keV. A 
subsequent paper will describe the soft X-ray spectrum of PG 1211+143 obtained with the 
LETG (Reeves et al. 2005, in preparation). 

In this paper we concentrate on the spectrum in the iron K band above 2 keV. A 
single power-law model provides an adequate fit to the continuum above 2 keV, with a 
photon index of F = 1.99 ± 0.06. The fit statistic obtained is x^/dof = 116.4/93, which is 
marginally acceptable (the rejection probabihty is 5x10"^). Figure 1 shows the LETG counts 
spectrum above 2 keV, the lower panel showing the residuals of the data compared to 
the continuum model. Each bin represents the FWHM resolution of the LETG. Two strong 
negative deviations that contribute significantly towards the total x^ are observed near 4.2 
and 4.9 keV, due to two absorption lines in the PG 1211+143 spectrum. To test this, two 
unresolved io = 10 eV) Gaussian absorption lines were added to the continuum model. A 
significant improvement in the spectral fit was obtained (x^/dof = 83.4/89), with addition of 
two absorption lines at 4.22 + 0.03 keV and 4.93 ±0.03 keV (observed frame) with equivalent 
widths of 35 ± 16 eV and 57 ± 23 eV. The improvement in fit statistic upon adding each fine 
was Ax^ = 13.8 and Ax^ = 19.9 respectively, while the overall fit statistic is now acceptable. 
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To constrain the width of the two absorption hnes, the LETG spectrum was rebinned to 
their maximum resolution, with each bin representing 1/4 FWHM resolution of the LETG 
(AA = 0.0125 A or Ai? = 16 eV at 4 keV). The rest-frame spectrum at this resolution 
near the absorption lines is shown in Figure 2. The spectral fit was minimized using the C- 
statistic (Cash 1979), due to the lower number of counts per bin at this resolution. The best 
fit width is cr ~ 21 eV (at 4.9 keV), corresponding to a FWHM velocity of ~ 3000 kms^^. 
However the fine width is poorly constrained, the 90% upper limit being < 7800 km s~^. 
The variabihty of the absorption hnes was tested by comparing the LETG spectra across 
all three Chandra orbits. The first LETG orbit was simultaneous with XMM-Newton, so 
the LETG and EPIC-pn spectra were fitted jointly. The results are shown in Table 1, while 
the C-statistic was used to minimize the spectral fits. The 4.2 and 4.9 kcV lines arc only 
formally detected during the 2nd orbit, although orbits 2 and 3 are statistically consistent. 
The lines do appear to vary in equivalent width between orbits 1 and 2 at > 90% confidence. 



3.1. Statistical Significance of the Absorption Lines 

In order to check that the absorption is not due to calibration effects, the PG 1211+143 
spectrum was compared to all the archived LETC observations of 3C273. No residuals 
are present in the summed or individual 3C273 spectra near the line energies 4.22 keV and 
4.93 kcV (Figure 3). The limit on the equivalent width of any absorption lines is < 4 eV in the 
summed spectrum. The source and background spectral extraction regions of PG 1211+143, 
in each of the three Chandra orbits, were also inspected. No anomalous pixels or rows were 
observed in the CCD array and no background sources were present within the extraction 
regions. The separate ±1 order spectra of PG 1211+143 were also found to be consistent. 
Thus the lines are not due to systematic or calibration effects in the LETG. 

In order to assess the statistical significance of the hnes, we wish to determine the 
chance probability of falsely detecting two lines, over an energy band of interest. In order to 
estimate this, 1000 LETG spectra were simulated, under the null hypothesis assumption that 
the spectrum consists of a broken power-law continuum from 0.4-8.0 keV, with no spectral 
lines. The distribution of Ax^ values obtained by fitting a single Gaussian line to each 
of the 1000 spectra simulated in the null hypothesis case can then be compared to the 
measured Ax^ in the real dataset to calculate the false probability of detection (Porquct ct 
al. 2004). The 1st order LETG spectra were simulated using identical continuum parameters 
and exposure as in the real PG 1211+143 dataset and binned at the FWHM resolution of the 
LETG. A single narrow Gaussian line (either in absorption or emission) was added to the 
continuum model and then fitted to each simulated spectrum. The line energy was stepped in 
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increments of one FWHM resolution bin over the 3-7 keV band and the simulated spectrum 
was re-fitted at each increment to minimize ^ . The minimum ^ obtained from adding the 
line to the spectral fit was compared to the obtained for the continuum model alone for 
each spectrum, in order to compute the distribution of Ax^ versus probability for all 1000 
simulated spectra. The 3-7 keV band was chosen for the line search as it is most likely that 
a strong iron Ka line will be observed over this energy range; i.e. < 3keV is outside the iron 
K band-pass for any feasible redshift, while the LETG has little effective area above 7 keV. 
The improvement in fit statistic obtained for adding the two lines to the actual dataset was 
Ax^ = 13.8 and Ax^ = 19.9 for each line. The simulations show that the false probability 
of detecting a line with equal or greater Ax^ is 0.034 and 0.004 for the 4.2 and 4.9 keV 
lines respectively. The probability of both lines being false detections is the product of these 
probabilities, which is 1.36 x 10~^. 

As an independent check, we calculate the line significances using the actual errors on 
the fitted parameters. The absorption line fluxes at 4.2 and 4.9 keV are —(4.9 ± 1.3) x 
10~^ photons cm~^ s"''' and — (5.9± 1.5) x 10^*^ photons cm~^ s"-*" respectively, while the con- 
tinuum fluxes at these energies are 1.40 ± 0.10 X 10"^ photons cm ^s ^ keV ^ and 1.04 ± 
0.08 X lO^'^ photons cm~^ s~^ keV~^. The equivalent widths of the lines are then —35 ± 10 eV 
and —57 ± 15 eV respectively. The errors are lex values and the line flux and continuum 
errors have been propogated. Thus the two absorption lines are significant at 3.5(7 and 3.8(j, 
with corresponding null probabihties of 4.6 x 10~^ and 1.4 x 10~^. However these are single 
trial probabilities for a line at an a-priori known energy. Taking into account the number of 
bins N searched for line features in the 3-7 keV range, then the null detection probabihty P 
will be P = 1 - (1 - Pi)^, where Pi is the single trial probability from above. As there are 
48 resolution bins in the 3-7 kcV range, the null probabilities for the 4.2 and 4.9 keV lines 
are 0.022 and 0.0067 respectively, in excellent agreement with the Monte-Carlo estimates. 



4. Discussion 

Two absorption fines are observed at 4.22 and 4.93 keV in the PG 1211+143 LETG 
spectrum, with a low probability that both are false detections. Although the lines are 
likely to be associated with redshifted Fe Ka, we first consider whether the lines could 
originate from other elements at lower velocities. The absorption line rest frame energies arc 
4.56 ± 0.03 keV and 5.33 ± 0.03 kcV. The closest transitions to those energies are Sc XXI Ka 
(at 4.53 keV) or V xxill (5.43 keV) respectively. The abundances of Sc or V are extremely 
low, ~ 1000 X less than Fe. Even if the abundances of were enhanced through spallation 
of iron nuclei by energetic protons on the disk surface (Skibo 1997), the strengths of these 
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lines would still be considerably weaker (<< 10%) than due to Fe Ka. One other possibihty 
is that the absorption originates from lower Z dements such as Si, S or Ca in a relativistic 
outflow. Blue-shifted lines with v ~ 0.1c have been detected in PG 1211+143 in a previous 
XMM-Newton observation by Pounds et al. (2003). In this case Ca XX at 4.10 keV would 
appear the most reasonable identification for the lower energy line, requiring an outfiow 
velocity of 32000 km s^^ However the Ca XX line is still likely to be weak, as the abundance 
ratio is Fe/Ca ~ 20 (Anders & Grevesse 1989). Si or S are more abundant than Ca, however 
if the fines originate from Si xiv Ka and S xvi Ka, then the outflow velocities are 0.68c and 
0.62c respectively. The kinetic power of such a fast flow would be a factor of ~ 40 greater 
than calculated by Pounds et al. (2003), probably exceeding 10^^ erg s~^ and the bolometric 
output of the quasar. Such a highly energetic flow seems implausible, unless the opening 
angle of the outflow is small. 

The most likely scenario is that the fines originate from highly redshifted iron Ka. 
Indeed there have been a small number of previous claims of redshifted Fe K absorption lines 
in other AGN. Nandra et al. (1999) first claimed tfie possibility of redshifted iron K-shell 
absorption, through ASCA observations of the Scyfcrt 1 NGC 3516. However recent Chandra 
and XMM-Newton observations (Turner et al. 2002) have showed that this is probably due 
to a strong narrow Ka core at 6.4 keV and variable emission lines at < 6.4 keV. Other 
claimed examples include a 6.2 keV fine in the quasar IE 1821-1-643 with Chandra HETG 
(Yaqoob & Serlemitsos 2005) and a weak redshifted fine in the quasar Q 0056-363 (Matt et 
al. 2005). Most recently the detection of transient redshifted Fe absorption fines in Mrk 509 
was claimed in Beppo-SAX data, with a velocity shift oi v — 0.21c (Dadina et al. 2005). 

If the lines in PG1211-M43 correspond to H-like iron (i.e. Fe xxvi Ka at 6.97 keV) 
then the redshifted velocities are 0.26c and 0.40c, while if the lines correspond to He-fike 
iron (Fe xxv Ka resonance at 6.70 keV), the velocities are 0.22c and 0.37c. Note that the 
separation of the two lines requires two distinct velocity components, as the higher order 
Fe K/3 lines are undetectable in this spectrum. It is unlikely that the lines result from a 
lower ionization than Fe XXV, as then a series of strong L-shell Fe lines and edges would be 
detected at lower energies, which is not the case. Furthermore the absorption lines cannot 
be from near neutral Fe (Fe l-xvil) at 6.4 kcV, because no strong Ka absorption line would 
be observed, as the L-shell would be fully populated. 

The two possible causes for the relativistic shifts are either gravitational redshift of 
photons near the massive black hole or infall of matter onto the black hole. In the former 
scenario, the absorbing matter would have to be located at 3.5 ± 0.5Rg and 4.8 ± 0.5i?g 
respectively from a Schwarzschild black hole (where i?g = GM/ is one gravitational radius) . 
This is within the last stable orbit around a Schwarzschild black hole and requires a Kerr 
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metric (Thorne 1974) for the matter to exist in a stable orbit within 6Rg. For a maximal pole- 
on Kerr black hole, the radii derived are then 3.2Rg and A.QRg respectively. Alternatively, 
the matter may be infalling onto the black hole, in which case a maximum velocity of ~ 0.38c 
(in the frame of the observer) can be achieved at 6Rg for a Schwarzschild metric. 

Modeling the two absorption line systems with Xstar (Kallman et al. 1996), requires a 
large column density of A^h = 4.0lfgX 10^^ cm~^ (assuming a turbulence of = 1000 km s^^ 
and solar abundances). The ionization parameter is log,^ = 3.9±0.3 erg cm s~^, with most of 
the absorption arising from Fc XXVI Ka. The equivalent column density of Fc XXVI is ~ 7 x 
10^*^ cm~^, while the optical depth at the K-shell edge is r < 0.1, i.e. the matter is optically- 
thin to the continuum. The velocity of the two absorption systems are v = (0.28 ± 0.01)c 
and V — (0.42 ± 0.01)c, while the line of sight covering fraction is > 50%. From the above 
parameters we can derive the physical properties of the gas. A black hole mass of 4 x 10^ Mq 
is assumed (Kaspi et al. 2000), an ionizing luminosity of L ~ 10^^ erg s"^ and a radial distance 
of ~ 10^^ cm (i.e. a few Rg). The ionization parameter is defined as ^ = L/nR^ (n is the 
electron density), while one can define a filling factor f as f = AR/R = N^/nR. A density 
of n ~ lO^^cm"^ is calculated, with a filling factor of ~1% and a cloud size of ~ 10^^ cm. 
Note that these should be considered order of magnitude estimates, due to the uncertainties 
in R and black hole mass. The matter may be clumped, perhaps from density perturbations 
on the disk or from localized filaments within the infalling matter. 

A fast outflow of ~ 0.1c has already been detected in an earlier XMM-Ncwton observa- 
tion of PG 1211 + 143 (Pounds et al. 2003) which is conflrmed from the analysis of the soft 
X-ray LETG data (Reeves et al. 2005, in preparation). One possibility is that part of this 
outflow does not escape the gravitational potential of the black hole (Murray & Chiang 1998; 
Proga, Stone, & Kallman 2000). If material is launched with v — 0.1c in PG 1211-1-143, then 
matter within a radius of it! < 100i?g may not escape the system. Thus it is possible to 
produce both red and blue-shifted hues in PG 1211-1-143, indeed in Mrk 509 both inflow and 
outflow also appear to be present (Dadina et al. 2005). From above, the mass infall rate 
required is M ~ ti fnR^vm^ ~ TiLvfm^/^ ~ lO^^gs"^ or ~ 0.1 Moyr"^. The kinetic power 
is then a few xlO^"^ erg s~^ for an infall velocity of f = c/3. Interestingly this is similar 
to the X-ray luminosity of PG 1211-1-143. Indeed a similar "aborted jet" model has been 
proposed by GhiseUini Haardt & Matt (2004), where matter in an outflow or jet fails to 
escape the gravitational potential well. The hard X-ray emission can be produced by shocks 
and collisions within the flow. 

Alternatively the absorption may originate from matter near the highly ionized inner 
disk surface, with the redshift being largely gravitational in origin. Ruszkowski & Fabian 
(2000) proposed such a model to explain the absorption line in NGC 3516, although the hues 
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will be broad unless the matter is clumped. Monitoring the line variability will help to reveal 
its origin, for instance it may possible to track the infall of a dense clump of matter towards 
the black hole. This bodes well for future longer observations of PG 1211+143, with XMM- 
Newton or Chandra, to track the line variability. Regardless of the mechanism, the large 
redshift and high covering fraction imply that most of the X-ray emission in PG 1211+143 
originates from a compact region within a few Rg of the black hole. 
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Fig. 1. — The upper panel shows the first order Chandra LETG spectrum of PG 1211+143 
above 2 keV; data points are shown as filled circles and la error bars as vertical lines. The 
best fitting continuum model, folded through the LETG response, is shown as a dashed line. 
The spectrum has been binned to the FWHM LETG spectral resolution and is plotted in the 
observed frame. The lower panel plots the residuals between the data and the continuum 
model, negative residuals correspond to a deficit of counts. Two absorption features are 
apparent at 4.2 and 4.9 keV respectively. 
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Fig. 2.— LETG spectrum of PG 1211+143 between 4-6 keV, binned at the full LETG 
resolution of AA = 0.0125 A (equivalent to 1/4 FWHM resolution per bin). Data are shown 
as filled circles, with dashed lines showing the la error bars and the best fit absorption 
line model as a solid line. The spectrum has been transformed into the rest frame of PG 
1211+143. If identified with Fe XXVI Ka, the redshift of each line would correspond to 
velocities of 0.26c and 0.40c. Note that the relative separation of the lines requires that two 
separate velocity components are present. 
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Fig. 3.— Data/Model ratio residuals of the PG 1211+143 (circles) and 3C 273 (stars) LETG 
spectra to power-law continuum models fitted above 2 keV. The data were binned at the 
FWHM resolution of the LETG. The absorption lines apparent between 4-5 keV in PG 
1211+143 do not appear in the 3C 273 spectrum, ruling out calibration effects. 
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Table 1. Table of iron line spectral fits. 



Parameter 


Mean 


Obsl 


Obs2 


Obs3 


Date^ 


21/06/2004 


21/06/2004 


23/06/2004 


25/06/2004 


Time^ 


03:09:04 


03:09:04 


03:57:43 


18:30:27 


Exposure^ 


133.6 


42.7 


48.2 


42.7 


pc 


1.99 ±0.06 


1.87 ±0.14 


2.00 ±0.10 


1.93 ±0.10 


2-10 


6.1 


4.0 


6.6 


7.3 


Energy'' 


4.22 ±0.03 


4.22^ 


4.22 ±0.03 


4.22^ 


EWs 


38 ±16 


< 25 


54 ±24 


< 46 


Energy'^ 


4.93 ±0.03 


4.93^ 


4.89 ± 0.04 


4.93^ 


EWs 


62 ±23 


< 21 


69 ±37 


< 57 



'^Start date and time (in UT) for observation. 
Exposure time in ks. 
'^Photon index above 2 keV. 
^2-10 keV flux in units lO"-*^^ ergs cm~^ s~^. 
*' Observed line energy in keV. 
^Denotes model parameter is fixed in fit. 
^Absorption line equivalent width in eV. 



